Relaxation of superfiow in a network: an application to the dislocation model of 

supersolidity of helium crystals 
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We have considered the dislocation network model for the supersolid state in 4 He crystals. In dif- 
ference with uniform two-dimensional and three-dimensional systems, the temperature of superfluid 
transition T c in the network is much smaller than the degeneracy temperature T<j. It is shown that a 
crossover into a quasi superfluid state occurs in the temperature interval between T c and Tij. Below 
the crossover temperature the time of decay of the flow increases exponentially under decrease in 
the temperature. The crossover has a continuous character and the crossover temperature does not 
depend on the density of dislocations. 

PACS numbers: 67.80.bd 



Forty years ago Andreev and Lifshitz 1 - predicted the 
possibility of realization of a ground state of quantum 
crystals with the number of sites larger than the number 
of atoms. Empty sites in such a state are called zero-point 
vacancies. In crystals consisting of Bose atoms zero-point 
vacancies behave as Bose quasiparticles and at low tem- 
peratures they become superfluid ones. Therefore the 
phenomenon predicted by Andreev and Lifshitz is called 
the supersolids. In 1970s the problem of supersolids was 
addressed a number of theoretical papers 2 -^^^, but nu- 
merous experimental efforts to discover this phenomenon 
failed (see Ref. 0) • The situation has changed completely 
after the experiment of Kim and Chan£ in which a non- 
classical rotational inertia (NCRI) of 4 He crystal was ob- 
served below certain critical temperature. NCRI means 
that the 4 He crystal does not rotate as a rigid body, but 
certain (superfluid) fraction of atoms is decoupled from 
the rotation. Further experiments^ ' 10 ' 11 ' 12 ' 13 confirmed 
the NCRI effect. It was established^ 1 ^ 1 ^ 1 ^ 4 . that the 
NCRI is not an intrinsic property of 4 He crystals. The 
amount of fraction decoupled from the oscillations (rota- 
tion) depends considerably on the degree of disorder in 
the crystal lattice, in particular, it becomes much smaller 
after annealing. A correlation between supersolid fea- 
tures and disorder was also observed in the direct flow 
experiment 1 ^. A common explanation of the correlation 
between disorder and supersolid behavior is that zero- 
point vacancies may emerge only at extended crystal lat- 
tice defects, presumable, in dislocation cores. The idea 
on superfluid behavior of dislocations was put forward in 
Ref. [l6| (originally, in a context of description of anoma- 
lous plastic properties of parahydrogen 17 ). In Ref. la 
this idea was applied for the explanation of NCRI. In[l8l 
and in the following 1 ^ papers superfluid properties of dis- 
locations were demonstrated by the first principles Monte 
Carlo calculations. A combined (dislocation network plus 
bulk) mechanism of supersolidity was considered in Ref. 
I20L The role of dislocations in this mechanism is simi- 
lar to one for the enhancement of superconductivity by 
dislocations^ 1 -. 

In this Rapid Communication we study the specifics 



of the superfluid transition in the network of one- 
dimensional (ID) wires and arrive at the following con- 
clusion. The temperature of the superfluid transition in 
the network is quite small and it depends on the length of 
segments of the network. At the same time, a crossover 
into a quasi superfluid state (a state with exponentially 
large time of relaxation of a flow) takes place at much 
larger temperatures, and the crossover temperature does 
not depend on the length of the segments. It is just 
the behavior observed in torsion experiments where the 
transition is continuous, and it occurs in the temperature 
interval independent of the NCRI fraction. 

The problem of the superfluid transition in a disloca- 
tion network was studied in Ref. |22j It was shown that 
in the two-dimensional (2D) network the critical tem- 
perature is proportional to the inverse length of the seg- 
ment. Let us derive the result of Ref. [22j using the fol- 
lowing simple arguments. For a 2D uniform medium the 
Berezinskii-Kosterlitz-Thouless (BKT) critical tempera- 
ture is determined by the equation T c = 7rh 2 n S 2(T c ) /2m, 
where n S 2 is the two-dimensional superfluid density. The 
BKT transition is connected with the energy of a vortex 
pair diverges logarithmically at large distances between 
the vortices. At such distances the difference between a 
uniform medium and a network is not important. There- 
fore, the critical temperature for the network is given by 
the BKT equation in which the uniform density n s % is 
replaced with the average 2D superfluid density n-2- The 
latter quantity can be expressed through the ID density 
of zero-point vacancies n and the length of the segment 
of the network I. For the quadratic network fii = 2n/l, 
and the critical temperature is T c — nh 2 n/ml (m is the 
effective mass of zero-point vacancies). 

In Ref. a microscopic model of Bose-Einstein con- 
densation (BEC) of non-interacting zero-point vacancies 
in a three-dimensional (3D) network was considered. It 
was shown that the BEC temperature depends signifi- 
cantly on the transparency of vertices (intersections of 
dislocations), and the highest value, reached at high 
transparency, is Tbec ~ 3h 2 n/2ml, so the superfluid crit- 
ical temperatures in 3D and 2D networks are of the same 
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order. 

The critical temperature T c in a network is much 
smaller than the degeneracy temperature for the ID Bosc 
gas of vacancies in a segment Td = h 2 n 2 /m (the small 
parameter is (nl) -1 -C 1)- In this respect the network 
differs from uniform 2D and 3D systems, where T c w Td- 
One can expect that the network may demonstrate quasi 
superfluid behavior in a wide range of temperatures be- 
tween T c and Td- 

To be more specific we consider a 2D network. In a 
uniform 2D system the vortex-antivortex pairs unbind 
above the BKT transition and vortices of opposite vor- 
ticities can move independently from each other. If a 
given vortex crosses the system in a direction perpendic- 
ular the flow, the superfluid phase difference along the 
flow changes on 2tt. In a uniform system a motion of 
a vortex across the flow is caused by the Magnus force 
and the viscous friction between the vortex and the nor- 
mal component. The network is a multiple connected 
system, the vortices correspond to circular currents, and 
they are pinned to given plaquettes. The vortex centers 
cannot move freely, but they can jump from one plaque- 
tte to another. This process becomes possible due to an 
emergence of phase slip (PS) centers at the segments. 

The theory of dissipation of supercurrent in ID chan- 
nels based on the idea of emergence of PS centers was 
put forward by Langer and Ambegaokar— . The main 
shortage of the theory 2 ^ is that it does not yield the 
pre-exponential factor in the expression for the relax- 
ation time. To obtain this factor one should consider 
the dynamics of transition of the system over the poten- 
tial barrier under appearance of a PS center. This prob- 
lem was solved in Ref. 25 on the base of the diffusive 
(time-dependent Ginzburg-Landau) equation. We solve 
the problem with the use of the wave (Gross-Pitaevskii) 
equation that describes an essentially different physical 
mechanism of relaxation of the supercurrent. 

If the vortex is already present in a given plaquette the 
PS of the proper sign provides annihilation of that vortex 
and creation of a vortex of the same vorticity in the neigh- 
bor plaquette. A vortex may jump to any neighbor pla- 
quette, but there is a preferable direction of such jumps 
in a system with a flow. Let us, for simplicity, consider 
a regular quadratic network with a flow directed parallel 
to the segments (say x direction). If the vortex is cen- 
tered in a given plaquette the superfluid velocities in the 
segments that form this plaquette read as va = v v + v s , 
vc = v v — v Sl vb = «d = v v , where A and C stand for the 
segments oriented along the x axis, while B and D - for 
the segments oriented along the y axis (perpendicular to 
the flow). Here v s is the flow velocity, and v v = nh/2ml 
is the contribution caused by the vortex. The frequency 
of PS in a segment is proportional to its length and it 
is a function of the superfluid velocity vi = lf(vi). One 
can see that vb = vd and in average the vortices do not 
move in the x direction (the direction of the flow). The 
PS frequencies for two other segments differ from each 
other and one extra jump in the perpendicular to the 



flow direction gains with the frequency 
Av = v A — v c = alv s , 



(1) 



where a = 2f'(v)\ v= o (here we imply the limit of small 
superfluid velocities). In the network of a rectangle shape 
of area S = L x x L y a vortex crosses the system with 
the frequency v C ross = N v Avl/L y , where N v = n v S is 
the total number of unbound vortices (n v is the vortex 
density). Each cross lowers the phase gradient on 2ir/L x 
that changes the flow velocity on Av s — —2irTi/mL x . The 
equation for v s , written in the differential form, reads as 



dv s 
dt 



,Av, = 



2nT%an v l 2 



(2) 



The solution of Eq. @ is v s — v^e^ 1 ! 1 ', where r = 
m / 2irfian v l 2 is the decay time. To compute a one 
should specify the mechanism of PS. We will describe 
the gas of zero-point vacancies in the dislocation core as a 
weakly non-ideal ID Bose gas with a complex condensate 
wave function (order parameter) ^(x,t) that satisfies the 
Gross-Pitaevskii (GP) equation. The ID GP equation 
with a repulsive point interaction has an exact solution 
that corresponds to a dark soliton. The dark soliton is a 
rarefaction that moves with a constant velocity u. 
The dark soliton is described by the function 26 



■ tanh 



ut 



(3) 



Here £ = h/mc is the coherence length, c = ^J^n/m 
is the sound velocity, 7 is the interaction constant, n — 
n(l — 2(£/£)yT — u 2 /c 2 ) -1 is the renormalized density 
(renormalization is the consequence of the conservation 
of the total number of zero-point vacancies) . The energy 
of the soliton reads as 



dx 



4 / u 
— Tine 1 — 



n 2 ) 

2\ 3/2 



(4) 



The soliton momentum p can be found by integration of 
equation dp = dEo/u: 



p = —2hn 




arcsm ■ 



C, 



(5) 



where C is the constant of integration. To determine C 
one can take into account that the soliton momentum 
is the difference of the momentum of the system with 
and without the soliton. The soliton emerges with the 
velocity u = +c — or u = — c + and at such u its 
momentum should be equal to zero. One can see that 
two conditions p u= ± c = yield two different integration 
constants C± = ±hmr. Therefore, one should consider 
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two species of the solitons (the "+" and "-" ones) with 
the momenta p± defined by Eq. ([5]) with C = C± . 

It is important to emphasize that these two species 
correspond two physically distinct solitons. According to 
Eq. the phase at the soliton has the additional shift 
A<p(u) — — 2 cot" 1 (it/Vc 2 — u 2 ). In any multiple con- 
nected system the phase satisfies the Onsager-Feynman 
quantization condition, and the appearance of a soliton 
should be accompanied by a change in the net velocity: 
v = vq + Av±(u). The function A<p(u) is discontinu- 
ous at u = with the jump equal 2tt. The functions 
Av±(u) should be continuous because small variation in 
u cannot result in a finite change in the net velocity. 
Since Av + (+c) = and Au_(— c) = 0, the function 
Av + (u) Av-(u). For instance, for a ID ring with 
the perimeter I we find Av± = h(±ir — 2 arcsin(u/c))/mZ 
(we take into account that Aip + mlAv±/h = mod 2ir). 
Thus, solitons of distinct species differ from each other 
by the change in the net velocity they induce. 

The soliton may change its velocity due to the inter- 
action with phonons and impurities. Since the soliton 
changes the net velocity in a continuous way, it is more 
consistent to consider PS as an entire process of creation 
of the soliton at u = +c and its annihilation at u = — c, 
or vise versa. The first possibility corresponds to the " +" 
solitons, and the second one - to the "-" ones. 

To obtain the frequency of PS we consider solitons as 
classical particles whose distribution functions f±(p,t) 
satisfy the Fokker- Planck equation 



yields 



df± 
dt 



ds± 
dp 



(0) 



where s± are the soliton fluxes in the momentum space. 
They read as 



d 

s± = A±f±- B± — f±. 

op 



(7) 



The coefficients A± and B± satisfy the relation A±fo t ± — 
B±dfo,±/dp = 0, where / ,± = ex p{—E±/T) is the equi- 
librium distribution function, and E± = Eg+p±v are the 
soliton energies at nonzero net velocity. Here we consider 
the case of slow relaxation. In this case one can neglect 
the explicit time dependence of f± and consider the fluxes 
s± as constant quantities. Using the relation between the 
coefficient A and B we rewrite Eq. ([7|) in the form 



o f d ( f ± 
dp V/o,± 



(8) 



The distribution of the solitons with small p is close to 
equilibrium one. The solitons with p± — > ±2irhn emerge 
only due to nonzero s±. To attain such a momentum 
the soliton should overcome the energy barrier AE = 
AhncfS. Therefore the fluxes are small and the distribu- 
tion functions f± at p = ±2irhn are much less than the 
equilibrium ones. The integration of © with the bound- 
ary conditions (/// )| P =o = 1 and (///o)|p=±2irRn = 



dp 



B±fo,± 



(9) 



Here we imply the case of small temperatures T <C AE 
and small net velocities v <C c. Then in the leading order 
the integral §§§ is evaluated as 

1 /2 

Bq ± / 2hnc\ ( Ahnc Trhnv±\ 



3T 



— 



(10) 

where i?o.± is the 'diffusion" coefficient B± at u = 0. Us- 
ing the exact form of fo one finds that B± = —A±T/ (u + 
v). The coefficient A = dp/dt is just the viscous friction 
force acting on solitons. At small velocities this force is 
proportional to the velocity of the soliton motion rela- 
tive the normal component A = —rj{u + v), that yields 
B = rjT, where rj is the friction coefficient (under as- 
sumption that rj is the same for "+" and "-" solitons, the 
coefficients A± = A and B± = B are the same as well). 

The soliton distribution functions are normalized by 
the condition n so i = (l/2Tth) f f dp, where n so i is the 
soliton density. The quantities s± are the fluxes in the 
direction of larger momenta p (negative sign of s_ means 
that actual direction of the flux is the opposite one) . The 
difference of their modules determines the frequency of 
PS's: 



-)l 



2nh 



(11) 



Using Eqs. (fTTj)) and (fTTjl one finds the parameter a and 
obtains the following expression for the decay time 



To- 



1 



T 



v l 2 \2nhnc 



1/2 



exp(^), (12) 



where t = 771/77. According to (TT2")) the crossover tem- 
perature is T = Ahnc/ 3. For weakly non- ideal Bose 
gas To/T d - y/-yn/T d < 1 and T /T c - l/£ > 1, so 
T c « T « T d . 

To evaluate the parameter To we consider the fric- 
tion connected with the interaction of the solitons with 
phonons. The dark soliton is the exact solution of the 
ID GP equation and it does not interact with phonons 
in ID. But the dislocation core is not a strict ID sys- 
tem. It is a quasi-lD system with a small, but finite 
cross-section. Such a system is described by an effec- 
tive ID GP equation with an additional higher order in 
* interaction term 2 !: iW^ /dt = -(h 2 /2m)d 2 ^/dx 2 + 
-y|\[r| 2 il/ _ ^l^l 4 ^'. Due to such a term the reflection co- 
efficient R ~ (71 77/7) 2 for the phonons that scatter on 
soliton is nonzero^. 

The parameter 7177 is evaluated as 7177 ~ 7(t'±/ / £) 2 , 
where r± is the radius the supefluid channel. As was 
shown in Ref. H3, at small u (and v — 0) the time deriva- 
tive of the soliton momentum is given by the expression 
p = —rju = —C(mRT/h)u, where the numerical factor 
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C ~ 1. It yields r ~ (h/T)(£,/r ± ) 4 . The scale of r is 
determined by the quantity h/T that is of order of 10~ 10 
c- 1 for T = 0.1 K. 

It is necessary to note one important point. One could 
think that since the GP equation we use is invariant with 
respect to translations it is impossible to describe the re- 
laxation of superflow in the GP approach. But our ap- 
proach is basically the same as commonly used for the 
computation of forces that act on a vortex in a 3D su- 
pcrfluid caused by its interaction with phonons and ro- 
tons (see, for instance Ref. 28). In the latter case the 
GP approach is applied for the computation of ampli- 
tudes of scattering of phonons on vortices and for finding 
the momentum flux over the cylindrical surface around 
the vortex line. In such a way one can obtain the rate 
of transfer of the momentum from a vortex to phonons, 
i.e.., to the normal component, which is assumed to be in 
equilibrium with the environment (walls, substrate, etc.) 

The specific of the ID system is that for the GP equa- 
tion with only cubic interaction term the phonons do 
not interact with solitons and there is no transfer of the 
momentum from solitons to the normal component. The 
fifth-order term switches on that interaction, and the mo- 
mentum transfers from solitons to ID phonons of the net- 
work. ID phonons interact with bulk phonons and due 
to such an interaction the normal component in the net- 
work remains in equilibrium with the crystal (since the 
relaxation time in the phonon subsystem is much smaller 
than t given by Eq. I|12p). Eventually, the momentum 
obtained from solitons is transferred to the crystal. Here 
we do not describe explicitly the mechanism of such a 
transfer, but just imply that the normal component is in 



equilibrium with the environment. 

Estimating (TT2")) we obtain that the decay time is of 
order of few seconds at T « 0.04Xb and it is of or- 
der of an hour at T w 0.03To. We note that at very 
small temperatures quantum jumps between states with 
different vortex configurations may become important 
(fo r ID rings the quantum jumps were studied in Refs. 
[29j|30h . We estimate the quantum correction to the r 
is To oc e n( - and this correction becomes important at 
T < = Tb/n£ -C To. The quantum correction re- 
sults only in a modification of the law of increase in the 
relaxation time under lowering of temperature and its ac- 
counting should not change the main conclusion on the 
emergence of the quasi superfluid state in the network 
well above the critical temperature. 

In this Rapid Communication we have considered quasi 
superfluidity in a 2D network. The situation in a 3D net- 
work should be qualitatively the same. In the latter case 
superflow may decay due to expansion or shrinking of 
vortex rings. The mechanism of expansion(shrinking) of 
vortex rings in a 3D network is basically the same as vor- 
tex motion in a 2D network: both of them are connected 
with the phase slips in segments of the network. 

In conclusion, we note that the results obtained can be 
also applied to multiple connected Bose-Einstein conden- 
sates of rarefied alkali gases in optical lattices^, where 
the quasi superfluid state can be observed directly. 
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